This is the analysis script to prepare the data to run analyses for the project “Anxiety, depression and trauma symptom change during the COVID-19 pandemic: retrospective versus objective assessment” - Young et al (2020)

Script written by K Purves, K Thompson, C Huebel and M Davies. Email: , , ,

#Set up

Clear global environment

Retrieve the current date to use it for file endings to not overwrite files when one exports files

#Packages Install packages (if they are not available in your version of R)

Load packages

Registered S3 method overwritten by 'pryr':
  method      from
  print.bytes Rcpp
For best results, restart R session and update pander using devtools:: or remotes::install_github('rapporter/pander')

Attaching package: 'polycor'
The following object is masked from 'package:psych':

    polyserial
corrplot 0.84 loaded
── Attaching packages ─────────────────────────────────────── tidyverse 1.3.0 ──
✓ ggplot2 3.3.2     ✓ purrr   0.3.4
✓ tibble  3.0.4     ✓ dplyr   1.0.2
✓ tidyr   1.1.2     ✓ stringr 1.4.0
✓ readr   1.3.1     ✓ forcats 0.5.0
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
x ggplot2::%+%()   masks psych::%+%()
x ggplot2::alpha() masks psych::alpha()
x dplyr::filter()  masks stats::filter()
x dplyr::lag()     masks stats::lag()
x tibble::view()   masks summarytools::view()

Colour palettes

Define colours for plotting this are the standard coping colours

Choose in this chunk which palette to use 04.07.2020 - Default to 2 colour COPING palette

ggplot theme

Set up ggplot theme for the plots

Data import

Source data file paths

Read in data

[1] 65223   643

Important vectors with collections of columns

Vectors for sum_score

 [1] "gad.sum_score_prepan"        "pcl.sum_score_prepan"       
 [3] "phq.sum_score_prepan"        "phq.sum_score_8items_prepan"
 [5] "ocir.sum_score_prepan"       "gad.sum_score_base"         
 [7] "gad.sum_score_retro"         "pcl.sum_score_base"         
 [9] "pcl.sum_score_retro"         "phq.sum_score_base"         
[11] "phq.sum_score_8items_base"   "phq.sum_score_retro"        
[13] "phq.sum_score_8items_retro"  "ocir.sum_score_base"        
[15] "ocir.sum_score_retro"       
[1] "phq.sum_score_prepan" "phq.sum_score_base"   "phq.sum_score_retro" 
[1] "gad.sum_score_prepan" "gad.sum_score_base"   "gad.sum_score_retro" 
[1] "ocir.sum_score_prepan" "ocir.sum_score_base"   "ocir.sum_score_retro" 
[1] "pcl.sum_score_prepan" "pcl.sum_score_base"   "pcl.sum_score_retro" 

Vectors for diff_score

 [1] "phq.diff_score_retro_prepan"  "phq.diff_score_base_prepan"  
 [3] "phq.diff_score_base_retro"    "gad.diff_score_retro_prepan" 
 [5] "gad.diff_score_base_prepan"   "gad.diff_score_base_retro"   
 [7] "pcl.diff_score_base_prepan"   "ocir.diff_score_retro_prepan"
 [9] "ocir.diff_score_base_prepan"  "ocir.diff_score_base_retro"  
[1] "phq.diff_score_retro_prepan" "phq.diff_score_base_prepan" 
[3] "phq.diff_score_base_retro"  
[1] "gad.diff_score_retro_prepan" "gad.diff_score_base_prepan" 
[3] "gad.diff_score_base_retro"  
[1] "ocir.diff_score_retro_prepan" "ocir.diff_score_base_prepan" 
[3] "ocir.diff_score_base_retro"  

Exclusions

Exclude GLAD/EDGI individuals that didn’t take part in COPING Exclude individuals with NA for the MHD categories Individuals with NAs for age/gender have already been excluded in pre-processing script, but they get excluded here to be safe and to double check previous programming.

[1] 65223   643
[1] 34465   643
[1] 30758

frequencies and descriptives in exclusion data

pandemic change gad

Frequencies  
dat$gad.pandemic_felt_feelings  
Type: Factor  

                         Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
--------------------- ------- --------- -------------- --------- --------------
           Much worse    5720     16.62          16.62     16.60          16.60
       A little worse   12774     37.12          53.75     37.06          53.66
         No different   13640     39.64          93.39     39.58          93.24
      A little better    1558      4.53          97.91      4.52          97.76
          Much better     718      2.09         100.00      2.08          99.84
                 <NA>      55                               0.16         100.00
                Total   34465    100.00         100.00    100.00         100.00

pandemic change phq

Frequencies  
dat$phq.pandemic_felt_feelings  
Type: Factor  

                         Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
--------------------- ------- --------- -------------- --------- --------------
           Much worse    5481     15.98          15.98     15.90          15.90
       A little worse   13745     40.07          56.04     39.88          55.78
         No different   12470     36.35          92.39     36.18          91.97
      A little better    1799      5.24          97.64      5.22          97.19
          Much better     811      2.36         100.00      2.35          99.54
                 <NA>     159                               0.46         100.00
                Total   34465    100.00         100.00    100.00         100.00

Ethnicity

Frequencies  
dat$Ethnicity_collapsed  
Type: Factor  

                                          Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
-------------------------------------- ------- --------- -------------- --------- --------------
                              European   30201     95.82          95.82     87.63          87.63
      Mixed or multiple ethnic origins     571      1.81          97.63      1.66          89.28
                Asian or Asian British     394      1.25          98.88      1.14          90.43
            African or African British     159      0.50          99.39      0.46          90.89
                                 Other     193      0.61         100.00      0.56          91.45
                                  <NA>    2947                               8.55         100.00
                                 Total   34465    100.00         100.00    100.00         100.00

Gender

Frequencies  
dat$Gender_collapsed  
Type: Factor  

                                 Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------------------------- ------- --------- -------------- --------- --------------
                         Male    9698     28.14          28.14     28.14          28.14
                       Female   24387     70.76          98.90     70.76          98.90
      Non-binary/Self-defined     380      1.10         100.00      1.10         100.00
                         <NA>       0                               0.00         100.00
                        Total   34465    100.00         100.00    100.00         100.00

age categorical

Frequencies  
dat$age_category_collapsed  
Type: Factor  

               Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------- ------- --------- -------------- --------- --------------
      16-18    1055      3.06           3.06      3.06           3.06
      19-25    2617      7.59          10.65      7.59          10.65
      26-35    5112     14.83          25.49     14.83          25.49
      36-45    4781     13.87          39.36     13.87          39.36
      46-55    6988     20.28          59.63     20.28          59.63
      56-65    8015     23.26          82.89     23.26          82.89
      66-70    3189      9.25          92.14      9.25          92.14
      71-75    2010      5.83          97.97      5.83          97.97
        76+     698      2.03         100.00      2.03         100.00
       <NA>       0                               0.00         100.00
      Total   34465    100.00         100.00    100.00         100.00

Employment

Frequencies  
dat$employment_prior_covid  
Type: Character  

                    Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
---------------- ------- --------- -------------- --------- --------------
        Employed    8923     27.14          27.14     25.89          25.89
      Key worker   12559     38.20          65.34     36.44          62.33
         Retired    8034     24.44          89.78     23.31          85.64
         Student    1507      4.58          94.36      4.37          90.01
      Unemployed    1854      5.64         100.00      5.38          95.39
            <NA>    1588                               4.61         100.00
           Total   34465    100.00         100.00    100.00         100.00

anxiety disorders

Frequencies  
dat$anxiety_disorders  
Type: Factor  

                             Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
------------------------- ------- --------- -------------- --------- --------------
      No anxiety disorder   17177     52.97          52.97     49.84          49.84
         Anxiety disorder   15251     47.03         100.00     44.25          94.09
                     <NA>    2037                               5.91         100.00
                    Total   34465    100.00         100.00    100.00         100.00

Depressive disorders

Frequencies  
dat$depressive_disorders  
Type: Factor  

                                Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
---------------------------- ------- --------- -------------- --------- --------------
      No depressive disorder   15261     45.77          45.77     44.28          44.28
         Depressive disorder   18084     54.23         100.00     52.47          96.75
                        <NA>    1120                               3.25         100.00
                       Total   34465    100.00         100.00    100.00         100.00

Depression and anxiety

Frequencies  
dat$depression_and_anxiety  
Type: Factor  

                                           Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
--------------------------------------- ------- --------- -------------- --------- --------------
      No depressive or anxiety disorder   12892     39.76          39.76     37.41          37.41
        Depressive and anxiety disorder   12963     39.97          79.73     37.61          75.02
               Only depressive disorder    4285     13.21          92.94     12.43          87.45
                  Only anxiety disorder    2288      7.06         100.00      6.64          94.09
                                   <NA>    2037                               5.91         100.00
                                  Total   34465    100.00         100.00    100.00         100.00

Eating disorders

Frequencies  
dat$eating_disorders  
Type: Factor  

                            Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
------------------------ ------- --------- -------------- --------- --------------
      No eating disorder   30748     93.54          93.54     89.22          89.22
         Eating disorder    2122      6.46         100.00      6.16          95.37
                    <NA>    1595                               4.63         100.00
                   Total   34465    100.00         100.00    100.00         100.00

OCD

Frequencies  
dat$obsessive_compulsive_disorders  
Type: Factor  

                                          Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
-------------------------------------- ------- --------- -------------- --------- --------------
      No obsessive compulsive disorder   30948     92.81          92.81     89.80          89.80
         Obsessive compulsive disorder    2397      7.19         100.00      6.95          96.75
                                  <NA>    1120                               3.25         100.00
                                 Total   34465    100.00         100.00    100.00         100.00

Psychosis

Frequencies  
dat$psychotic_disorders  
Type: Factor  

                               Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
--------------------------- ------- --------- -------------- --------- --------------
      No psychotic disorder   32239     98.20          98.20     93.54          93.54
         Psychotic disorder     590      1.80         100.00      1.71          95.25
                       <NA>    1636                               4.75         100.00
                      Total   34465    100.00         100.00    100.00         100.00

Bipolar

Frequencies  
dat$mhd.mania_hypomania_bipolar_or_manicdepression  
Type: Factor  

                                                           Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
------------------------------------------------------- ------- --------- -------------- --------- --------------
      Not Mania, hypomania, bipolar or manic-depression   32156     96.43          96.43     93.30          93.30
          Mania, hypomania, bipolar or manic-depression    1189      3.57         100.00      3.45          96.75
                                                   <NA>    1120                               3.25         100.00
                                                  Total   34465    100.00         100.00    100.00         100.00

PTSD

Frequencies  
dat$mhd.posttraumatic_stress_disorder_ptsd  
Type: Factor  

                                                   Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------------------------------------------- ------- --------- -------------- --------- --------------
      Not Post-traumatic stress disorder (PTSD)   30556     91.64          91.64     88.66          88.66
          Post-traumatic stress disorder (PTSD)    2789      8.36         100.00      8.09          96.75
                                           <NA>    1120                               3.25         100.00
                                          Total   34465    100.00         100.00    100.00         100.00

ASD

Frequencies  
dat$autism_spectrum_disorder  
Type: Factor  

                                     Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
--------------------------------- ------- --------- -------------- --------- --------------
      No autism spectrum disorder   32202     98.09          98.09     93.43          93.43
         Autism spectrum disorder     627      1.91         100.00      1.82          95.25
                             <NA>    1636                               4.75         100.00
                            Total   34465    100.00         100.00    100.00         100.00

ADHD

Frequencies  
dat$mhd.attention_deficit_hyperactivity_disorder  
Type: Factor  

                                                                                          Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
-------------------------------------------------------------------------------------- ------- --------- -------------- --------- --------------
      Not Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD)   32480     98.94          98.94     94.24          94.24
          Attention deficit or attention deficit and hyperactivity disorder (ADD/ADHD)     349      1.06         100.00      1.01          95.25
                                                                                  <NA>    1636                               4.75         100.00
                                                                                 Total   34465    100.00         100.00    100.00         100.00

Personality disorders

Frequencies  
dat$mhd.personality_disorder  
Type: Factor  

                                  Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
------------------------------ ------- --------- -------------- --------- --------------
      Not Personality disorder   31501     95.95          95.95     91.40          91.40
          Personality disorder    1328      4.05         100.00      3.85          95.25
                          <NA>    1636                               4.75         100.00
                         Total   34465    100.00         100.00    100.00         100.00

Controls

Frequencies  
dat$control  
Type: Factor  

                    Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
---------------- ------- --------- -------------- --------- --------------
      No control   21348     64.16          64.16     61.94          61.94
         Control   11925     35.84         100.00     34.60          96.54
            <NA>    1192                               3.46         100.00
           Total   34465    100.00         100.00    100.00         100.00

Dont know or prefer not to say

Frequencies  
dat$DidNotRespondMHD  
Type: Factor  

                           Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------------------- ------- --------- -------------- --------- --------------
      Prefer not to say    1086      3.15           3.15      3.15           3.15
       Responded to MHD   33379     96.85         100.00     96.85         100.00
                   <NA>       0                               0.00         100.00
                  Total   34465    100.00         100.00    100.00         100.00

Data availability per questionnaire GAD, PHQ, PCL, and OCI-R

++MD: These ‘data availability’ chunks caused my R to crash, so I’m hashing them out for now. Others can test if it works for them and unhash as needed, otherwise I think they should be removed.

Data availability GAD

Data availability PHQ

Data availability PCL

Data availability OCI-R

Whole sample

Participants per Sample

Eating Disorders Genetics Initative Genetic Links to Anxiety and Depression National BioResource Rapid Assessment of MP

Frequencies  
dat$Sample  
Type: Character  

               Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------- ------- --------- -------------- --------- --------------
       EDGI      65      0.19           0.19      0.19           0.19
       GLAD   12653     36.71          36.90     36.71          36.90
        NBR   13096     38.00          74.90     38.00          74.90
       RAMP    8651     25.10         100.00     25.10         100.00
       <NA>       0                               0.00         100.00
      Total   34465    100.00         100.00    100.00         100.00

Psychiatric disorders

Frequencies  
dat$disorders_total_count  
Type: Numeric  

               Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------- ------- --------- -------------- --------- --------------
          0   12386     38.85          38.85     35.94          35.94
          1    5527     17.33          56.18     16.04          51.97
          2    8340     26.16          82.34     24.20          76.17
          3    3337     10.47          92.81      9.68          85.86
          4    1393      4.37          97.17      4.04          89.90
          5     569      1.78          98.96      1.65          91.55
          6     239      0.75          99.71      0.69          92.24
          7      64      0.20          99.91      0.19          92.43
          8      23      0.07          99.98      0.07          92.49
          9       6      0.02         100.00      0.02          92.51
       <NA>    2581                               7.49         100.00
      Total   34465    100.00         100.00    100.00         100.00

#HR

#SB - Impute no change prepandemic scores

 [1] 21 11 16 10  1  4 13  0 14 11  0  9 16  8  7 14 14  6  4  5  9  9 21 19  4
[26]  6  7  9  7 14
 [1] NA  2  9  5 NA 21 NA NA NA  3 NA  3 NA  3  3 13 NA NA  3  0 NA  5 15 17 NA
[26] NA  2 NA  4  4
 [1] No different   Much worse     A little worse A little worse No different  
 [6] Much better    No different   No different   No different   A little worse
5 Levels: Much worse A little worse No different ... Much better
 [1] 3 1 2 2 3 5 3 3 3 2
 [1] 21 11 16 10  1  4 13  0 14 11  0  9 16  8  7 14 14  6  4  5  9  9 21 19  4
[26]  6  7  9  7 14
 [1] 21  2  9  5  1 21 13  0 14  3  0  3 16  3  3 13 14  6  3  0  9  5 15 17  4
[26]  6  2  9  4  4
 [1] 22  9 15 11  0  9 23  5 23 11  1 20 18 15  4 11 15 19 12 10 17 11 24 27  6
[26]  9  3 11  3  8
 [1] 23  1  7  7 NA 20 18 NA 23  5 NA  3 NA  7  3 11 NA NA 12  2 NA  6 19 27 NA
[26] NA  1 NA NA  3
 [1] A little worse  Much worse      A little worse  A little worse 
 [5] No different    A little better A little worse  No different   
 [9] A little worse  A little worse 
5 Levels: Much worse A little worse No different ... Much better
 [1] 2 1 2 2 3 4 2 3 2 2
 [1] 22  9 15 11  0  9 23  5 23 11  1 20 18 15  4 11 15 19 12 10 17 11 24 27  6
[26]  9  3 11  3  8
 [1] 23  1  7  7  0 20 18  5 23  5  1  3 18  7  3 11 15 19 12  2 17  6 19 27  6
[26]  9  1 11  3  3

#Baseline survey completion (N per month)

 [1] "2020-04-30 20:51:40 UTC" "2020-04-30 18:55:12 UTC"
 [3] "2020-04-30 19:02:04 UTC" "2020-04-30 18:59:37 UTC"
 [5] "2020-05-01 10:02:34 UTC" "2020-04-30 22:16:01 UTC"
 [7] "2020-04-30 20:36:37 UTC" "2020-04-30 19:02:03 UTC"
 [9] "2020-04-30 18:57:33 UTC" "2020-05-01 08:54:01 UTC"
[11] "2020-04-30 21:47:03 UTC" "2020-04-30 19:10:42 UTC"
[13] "2020-04-30 19:29:36 UTC" "2020-04-30 19:20:42 UTC"
[15] "2020-05-01 15:09:23 UTC" "2020-04-30 18:59:55 UTC"
[17] "2020-05-04 16:40:21 UTC" "2020-05-01 16:06:45 UTC"
[19] "2020-05-19 10:34:16 UTC" "2020-04-30 19:11:26 UTC"

Vector

Correlation matrix

Warning in polychor(x, y, ML = ML, std.err = std.err): inadmissible correlation
set to -0.9999

Warning in polychor(x, y, ML = ML, std.err = std.err): inadmissible correlation
set to -0.9999
Warning in hetcor.data.frame(as.data.frame(dat[, mhd.items]), use =
"pairwise.complete.obs"): the correlation matrix has been adjusted to make it
positive-definite

Two-Step Estimates

Correlations/Type of Correlation:
                                               depressive_disorders
depressive_disorders                                              1
anxiety_disorders                                            0.8377
eating_disorders                                             0.5343
obsessive_compulsive_disorders                               0.5269
psychotic_disorders                                          0.5033
mhd.mania_hypomania_bipolar_or_manicdepression               0.5367
mhd.posttraumatic_stress_disorder_ptsd                       0.5807
autism_spectrum_disorder                                     0.3882
mhd.attention_deficit_hyperactivity_disorder                 0.3696
mhd.personality_disorder                                     0.6383
control                                                     -0.8096
                                               anxiety_disorders
depressive_disorders                                  Polychoric
anxiety_disorders                                              1
eating_disorders                                          0.4969
obsessive_compulsive_disorders                            0.6313
psychotic_disorders                                       0.4568
mhd.mania_hypomania_bipolar_or_manicdepression             0.431
mhd.posttraumatic_stress_disorder_ptsd                    0.5833
autism_spectrum_disorder                                  0.4385
mhd.attention_deficit_hyperactivity_disorder              0.4147
mhd.personality_disorder                                  0.5838
control                                                  -0.8075
                                               eating_disorders
depressive_disorders                                 Polychoric
anxiety_disorders                                    Polychoric
eating_disorders                                              1
obsessive_compulsive_disorders                           0.5707
psychotic_disorders                                      0.3641
mhd.mania_hypomania_bipolar_or_manicdepression           0.3653
mhd.posttraumatic_stress_disorder_ptsd                   0.3867
autism_spectrum_disorder                                 0.2962
mhd.attention_deficit_hyperactivity_disorder             0.3309
mhd.personality_disorder                                 0.5191
control                                                 -0.6858
                                               obsessive_compulsive_disorders
depressive_disorders                                               Polychoric
anxiety_disorders                                                  Polychoric
eating_disorders                                                   Polychoric
obsessive_compulsive_disorders                                              1
psychotic_disorders                                                    0.3917
mhd.mania_hypomania_bipolar_or_manicdepression                         0.3917
mhd.posttraumatic_stress_disorder_ptsd                                 0.4481
autism_spectrum_disorder                                               0.3976
mhd.attention_deficit_hyperactivity_disorder                           0.3908
mhd.personality_disorder                                                0.502
control                                                               -0.7346
                                               psychotic_disorders
depressive_disorders                                    Polychoric
anxiety_disorders                                       Polychoric
eating_disorders                                        Polychoric
obsessive_compulsive_disorders                          Polychoric
psychotic_disorders                                              1
mhd.mania_hypomania_bipolar_or_manicdepression              0.7086
mhd.posttraumatic_stress_disorder_ptsd                      0.4476
autism_spectrum_disorder                                    0.2901
mhd.attention_deficit_hyperactivity_disorder                0.3375
mhd.personality_disorder                                    0.6044
control                                                     -0.701
                                               mhd.mania_hypomania_bipolar_or_manicdepression
depressive_disorders                                                               Polychoric
anxiety_disorders                                                                  Polychoric
eating_disorders                                                                   Polychoric
obsessive_compulsive_disorders                                                     Polychoric
psychotic_disorders                                                                Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression                                              1
mhd.posttraumatic_stress_disorder_ptsd                                                 0.3492
autism_spectrum_disorder                                                               0.2813
mhd.attention_deficit_hyperactivity_disorder                                            0.342
mhd.personality_disorder                                                               0.5855
control                                                                               -0.6882
                                               mhd.posttraumatic_stress_disorder_ptsd
depressive_disorders                                                       Polychoric
anxiety_disorders                                                          Polychoric
eating_disorders                                                           Polychoric
obsessive_compulsive_disorders                                             Polychoric
psychotic_disorders                                                        Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression                             Polychoric
mhd.posttraumatic_stress_disorder_ptsd                                              1
autism_spectrum_disorder                                                       0.2916
mhd.attention_deficit_hyperactivity_disorder                                   0.2851
mhd.personality_disorder                                                       0.5077
control                                                                       -0.6867
                                               autism_spectrum_disorder
depressive_disorders                                         Polychoric
anxiety_disorders                                            Polychoric
eating_disorders                                             Polychoric
obsessive_compulsive_disorders                               Polychoric
psychotic_disorders                                          Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression               Polychoric
mhd.posttraumatic_stress_disorder_ptsd                       Polychoric
autism_spectrum_disorder                                              1
mhd.attention_deficit_hyperactivity_disorder                     0.5726
mhd.personality_disorder                                         0.3777
control                                                         -0.6313
                                               mhd.attention_deficit_hyperactivity_disorder
depressive_disorders                                                             Polychoric
anxiety_disorders                                                                Polychoric
eating_disorders                                                                 Polychoric
obsessive_compulsive_disorders                                                   Polychoric
psychotic_disorders                                                              Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression                                   Polychoric
mhd.posttraumatic_stress_disorder_ptsd                                           Polychoric
autism_spectrum_disorder                                                         Polychoric
mhd.attention_deficit_hyperactivity_disorder                                              1
mhd.personality_disorder                                                             0.3698
control                                                                             -0.6354
                                               mhd.personality_disorder
depressive_disorders                                         Polychoric
anxiety_disorders                                            Polychoric
eating_disorders                                             Polychoric
obsessive_compulsive_disorders                               Polychoric
psychotic_disorders                                          Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression               Polychoric
mhd.posttraumatic_stress_disorder_ptsd                       Polychoric
autism_spectrum_disorder                                     Polychoric
mhd.attention_deficit_hyperactivity_disorder                 Polychoric
mhd.personality_disorder                                              1
control                                                         -0.7725
                                                  control
depressive_disorders                           Polychoric
anxiety_disorders                              Polychoric
eating_disorders                               Polychoric
obsessive_compulsive_disorders                 Polychoric
psychotic_disorders                            Polychoric
mhd.mania_hypomania_bipolar_or_manicdepression Polychoric
mhd.posttraumatic_stress_disorder_ptsd         Polychoric
autism_spectrum_disorder                       Polychoric
mhd.attention_deficit_hyperactivity_disorder   Polychoric
mhd.personality_disorder                       Polychoric
control                                                 1

Standard Errors/Numbers of Observations:
                                               depressive_disorders
depressive_disorders                                          33345
anxiety_disorders                                          0.004037
eating_disorders                                            0.01226
obsessive_compulsive_disorders                               0.0117
psychotic_disorders                                         0.02142
mhd.mania_hypomania_bipolar_or_manicdepression              0.01569
mhd.posttraumatic_stress_disorder_ptsd                      0.01078
autism_spectrum_disorder                                    0.02049
mhd.attention_deficit_hyperactivity_disorder                0.02606
mhd.personality_disorder                                    0.01498
control                                                         Inf
                                               anxiety_disorders
depressive_disorders                                       32428
anxiety_disorders                                          32428
eating_disorders                                         0.01209
obsessive_compulsive_disorders                           0.01054
psychotic_disorders                                       0.0202
mhd.mania_hypomania_bipolar_or_manicdepression           0.01546
mhd.posttraumatic_stress_disorder_ptsd                   0.01032
autism_spectrum_disorder                                 0.01985
mhd.attention_deficit_hyperactivity_disorder             0.02525
mhd.personality_disorder                                 0.01397
control                                                      Inf
                                               eating_disorders
depressive_disorders                                      32836
anxiety_disorders                                         31919
eating_disorders                                          32870
obsessive_compulsive_disorders                          0.01299
psychotic_disorders                                     0.02372
mhd.mania_hypomania_bipolar_or_manicdepression          0.01942
mhd.posttraumatic_stress_disorder_ptsd                  0.01538
autism_spectrum_disorder                                0.02513
mhd.attention_deficit_hyperactivity_disorder            0.02924
mhd.personality_disorder                                 0.0157
control                                                   1.851
                                               obsessive_compulsive_disorders
depressive_disorders                                                    33345
anxiety_disorders                                                       32428
eating_disorders                                                        32836
obsessive_compulsive_disorders                                          33345
psychotic_disorders                                                   0.02248
mhd.mania_hypomania_bipolar_or_manicdepression                        0.01825
mhd.posttraumatic_stress_disorder_ptsd                                0.01397
autism_spectrum_disorder                                              0.02216
mhd.attention_deficit_hyperactivity_disorder                           0.0268
mhd.personality_disorder                                              0.01562
control                                                                 11.73
                                               psychotic_disorders
depressive_disorders                                         32804
anxiety_disorders                                            31887
eating_disorders                                             32826
obsessive_compulsive_disorders                               32804
psychotic_disorders                                          32829
mhd.mania_hypomania_bipolar_or_manicdepression             0.01621
mhd.posttraumatic_stress_disorder_ptsd                     0.02076
autism_spectrum_disorder                                   0.03517
mhd.attention_deficit_hyperactivity_disorder               0.03917
mhd.personality_disorder                                   0.01887
control                                                      3.285
                                               mhd.mania_hypomania_bipolar_or_manicdepression
depressive_disorders                                                                    33345
anxiety_disorders                                                                       32428
eating_disorders                                                                        32836
obsessive_compulsive_disorders                                                          33345
psychotic_disorders                                                                     32804
mhd.mania_hypomania_bipolar_or_manicdepression                                          33345
mhd.posttraumatic_stress_disorder_ptsd                                                0.01847
autism_spectrum_disorder                                                              0.02949
mhd.attention_deficit_hyperactivity_disorder                                          0.03287
mhd.personality_disorder                                                              0.01632
control                                                                                 1.964
                                               mhd.posttraumatic_stress_disorder_ptsd
depressive_disorders                                                            33345
anxiety_disorders                                                               32428
eating_disorders                                                                32836
obsessive_compulsive_disorders                                                  33345
psychotic_disorders                                                             32804
mhd.mania_hypomania_bipolar_or_manicdepression                                  33345
mhd.posttraumatic_stress_disorder_ptsd                                          33345
autism_spectrum_disorder                                                      0.02388
mhd.attention_deficit_hyperactivity_disorder                                  0.02903
mhd.personality_disorder                                                      0.01514
control                                                                         1.983
                                               autism_spectrum_disorder
depressive_disorders                                              32804
anxiety_disorders                                                 31887
eating_disorders                                                  32826
obsessive_compulsive_disorders                                    32804
psychotic_disorders                                               32829
mhd.mania_hypomania_bipolar_or_manicdepression                    32804
mhd.posttraumatic_stress_disorder_ptsd                            32804
autism_spectrum_disorder                                          32829
mhd.attention_deficit_hyperactivity_disorder                    0.02798
mhd.personality_disorder                                        0.02533
control                                                           3.353
                                               mhd.attention_deficit_hyperactivity_disorder
depressive_disorders                                                                  32804
anxiety_disorders                                                                     31887
eating_disorders                                                                      32826
obsessive_compulsive_disorders                                                        32804
psychotic_disorders                                                                   32829
mhd.mania_hypomania_bipolar_or_manicdepression                                        32804
mhd.posttraumatic_stress_disorder_ptsd                                                32804
autism_spectrum_disorder                                                              32829
mhd.attention_deficit_hyperactivity_disorder                                          32829
mhd.personality_disorder                                                            0.03046
control                                                                               2.038
                                               mhd.personality_disorder control
depressive_disorders                                              32804   33271
anxiety_disorders                                                 31887   32382
eating_disorders                                                  32826   32800
obsessive_compulsive_disorders                                    32804   33271
psychotic_disorders                                               32829   32777
mhd.mania_hypomania_bipolar_or_manicdepression                    32804   33271
mhd.posttraumatic_stress_disorder_ptsd                            32804   33271
autism_spectrum_disorder                                          32829   32777
mhd.attention_deficit_hyperactivity_disorder                      32829   32777
mhd.personality_disorder                                          32829   32777
control                                                           1.749   33273

Create the correlation matrix plot using the corrplot package

#Correlations of retro impute and prepan ##Create function to use for all of the corplots

###Whole sample prepan retro correlation

Frequencies  
prepan_retro_dat$Sample  
Type: Character  

               Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------- ------- --------- -------------- --------- --------------
       EDGI      65      0.51           0.51      0.51           0.51
       GLAD   12653     99.49         100.00     99.49         100.00
       <NA>       0                               0.00         100.00
      Total   12718    100.00         100.00    100.00         100.00

Two-Step Estimates

Correlations/Type of Correlation:
                     gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan                    1             Pearson
gad.sum_score_retro                0.5143                   1
phq.sum_score_prepan               0.6904              0.4485
phq.sum_score_retro                0.4443              0.7445
                     phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan              Pearson             Pearson
gad.sum_score_retro               Pearson             Pearson
phq.sum_score_prepan                    1             Pearson
phq.sum_score_retro                0.5886                   1

Standard Errors/Numbers of Observations:
                     gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan                12692                9063
gad.sum_score_retro              0.007727                9082
phq.sum_score_prepan             0.004652            0.008394
phq.sum_score_retro              0.008274            0.004961
                     phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan                12662                9411
gad.sum_score_retro                  9058                8074
phq.sum_score_prepan                12681                9405
phq.sum_score_retro              0.006739                9431

P-values for Tests of Bivariate Normality:
                     gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan                                         
gad.sum_score_retro            6.754e-309                    
phq.sum_score_prepan           2.342e-223          2.296e-275
phq.sum_score_retro            4.022e-131          7.718e-221
                     phq.sum_score_prepan
gad.sum_score_prepan                     
gad.sum_score_retro                      
phq.sum_score_prepan                     
phq.sum_score_retro            1.101e-125
Warning in corCi(x = x, keys = keys, n.iter = n.iter, p = p, overlap = overlap, : 
Finding confidence intervals from a correlation matrix, but  n is not specified, arbitrarily set to 1000

Frequencies  
whole_retro_prepan_scatter$Sample  
Type: Character  

               Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------- ------- --------- -------------- --------- --------------
       EDGI      65      0.51           0.51      0.51           0.51
       GLAD   12653     99.49         100.00     99.49         100.00
       <NA>       0                               0.00         100.00
      Total   12718    100.00         100.00    100.00         100.00
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 33 rows containing non-finite values (stat_smooth).
Warning: Removed 33 rows containing missing values (geom_point).

Warning: Removed 33 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 42 rows containing non-finite values (stat_smooth).
Warning: Removed 42 rows containing missing values (geom_point).

Warning: Removed 42 rows containing missing values (geom_point).

###Little/much worse prepan retro correlation

[1] "gad.sum_score_prepan" "gad.sum_score_retro"  "phq.sum_score_prepan"
[4] "phq.sum_score_retro" 
Frequencies  
prepan_retro_dat_worse$Sample  
Type: Character  

              Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------- ------ --------- -------------- --------- --------------
       EDGI     42      0.54           0.54      0.54           0.54
       GLAD   7795     99.46         100.00     99.46         100.00
       <NA>      0                               0.00         100.00
      Total   7837    100.00         100.00    100.00         100.00

Two-Step Estimates

Correlations/Type of Correlation:
                     gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan                    1             Pearson
gad.sum_score_retro                0.5476                   1
phq.sum_score_prepan                 0.69              0.4785
phq.sum_score_retro                0.4592              0.7366
                     phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan              Pearson             Pearson
gad.sum_score_retro               Pearson             Pearson
phq.sum_score_prepan                    1             Pearson
phq.sum_score_retro                0.5976                   1

Standard Errors/Numbers of Observations:
                     gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan                 7820                7808
gad.sum_score_retro              0.007924                7825
phq.sum_score_prepan             0.005931            0.008729
phq.sum_score_retro              0.009446            0.005473
                     phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan                 7805                6982
gad.sum_score_retro                  7805                6988
phq.sum_score_prepan                 7817                6979
phq.sum_score_retro              0.007697                6997

P-values for Tests of Bivariate Normality:
                     gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan                                         
gad.sum_score_retro            8.019e-142                    
phq.sum_score_prepan           1.022e-144          2.245e-111
phq.sum_score_retro            1.966e-123          1.629e-103
                     phq.sum_score_prepan
gad.sum_score_prepan                     
gad.sum_score_retro                      
phq.sum_score_prepan                     
phq.sum_score_retro             2.61e-134
Warning in corCi(x = x, keys = keys, n.iter = n.iter, p = p, overlap = overlap, : 
Finding confidence intervals from a correlation matrix, but  n is not specified, arbitrarily set to 1000

Frequencies  
little_much_worse_retro_prepan$Sample  
Type: Character  

               Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------- ------- --------- -------------- --------- --------------
       EDGI      65      0.51           0.51      0.51           0.51
       GLAD   12653     99.49         100.00     99.49         100.00
       <NA>       0                               0.00         100.00
      Total   12718    100.00         100.00    100.00         100.00
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 14 rows containing non-finite values (stat_smooth).
Warning: Removed 14 rows containing missing values (geom_point).

Warning: Removed 14 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 15 rows containing non-finite values (stat_smooth).
Warning: Removed 15 rows containing missing values (geom_point).

Warning: Removed 15 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 29 rows containing non-finite values (stat_smooth).
Warning: Removed 29 rows containing missing values (geom_point).

Warning: Removed 29 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 11 rows containing non-finite values (stat_smooth).
Warning: Removed 11 rows containing missing values (geom_point).

Warning: Removed 11 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 21 rows containing non-finite values (stat_smooth).
Warning: Removed 21 rows containing missing values (geom_point).

Warning: Removed 21 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 32 rows containing non-finite values (stat_smooth).
Warning: Removed 32 rows containing missing values (geom_point).

Warning: Removed 32 rows containing missing values (geom_point).

###Little/much better prepan retro correlation

[1] "gad.sum_score_prepan" "gad.sum_score_retro"  "phq.sum_score_prepan"
[4] "phq.sum_score_retro" 
Frequencies  
prepan_retro_dat_better$Sample  
Type: Character  

              Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------- ------ --------- -------------- --------- --------------
       EDGI      7      0.56           0.56      0.56           0.56
       GLAD   1252     99.44         100.00     99.44         100.00
       <NA>      0                               0.00         100.00
      Total   1259    100.00         100.00    100.00         100.00

Two-Step Estimates

Correlations/Type of Correlation:
                     gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan                    1             Pearson
gad.sum_score_retro                0.4858                   1
phq.sum_score_prepan               0.6512              0.4084
phq.sum_score_retro                0.3567              0.6977
                     phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan              Pearson             Pearson
gad.sum_score_retro               Pearson             Pearson
phq.sum_score_prepan                    1             Pearson
phq.sum_score_retro                0.5346                   1

Standard Errors/Numbers of Observations:
                     gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan                 1257                1255
gad.sum_score_retro               0.02158                1257
phq.sum_score_prepan              0.01628             0.02355
phq.sum_score_retro               0.02651             0.01559
                     phq.sum_score_prepan phq.sum_score_retro
gad.sum_score_prepan                 1253                1085
gad.sum_score_retro                  1253                1086
phq.sum_score_prepan                 1255                1084
phq.sum_score_retro               0.02171                1087

P-values for Tests of Bivariate Normality:
                     gad.sum_score_prepan gad.sum_score_retro
gad.sum_score_prepan                                         
gad.sum_score_retro             2.716e-29                    
phq.sum_score_prepan            7.975e-34            3.25e-08
phq.sum_score_retro             1.024e-27           1.832e-07
                     phq.sum_score_prepan
gad.sum_score_prepan                     
gad.sum_score_retro                      
phq.sum_score_prepan                     
phq.sum_score_retro             1.275e-12
Warning in corCi(x = x, keys = keys, n.iter = n.iter, p = p, overlap = overlap, : 
Finding confidence intervals from a correlation matrix, but  n is not specified, arbitrarily set to 1000

Frequencies  
little_much_better_retro_prepan$Sample  
Type: Character  

               Freq   % Valid   % Valid Cum.   % Total   % Total Cum.
----------- ------- --------- -------------- --------- --------------
       EDGI      65      0.51           0.51      0.51           0.51
       GLAD   12653     99.49         100.00     99.49         100.00
       <NA>       0                               0.00         100.00
      Total   12718    100.00         100.00    100.00         100.00
Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 3 rows containing non-finite values (stat_smooth).
Warning: Removed 3 rows containing missing values (geom_point).

Warning: Removed 3 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 1 rows containing non-finite values (stat_smooth).
Warning: Removed 1 rows containing missing values (geom_point).

Warning: Removed 1 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 4 rows containing non-finite values (stat_smooth).
Warning: Removed 4 rows containing missing values (geom_point).

Warning: Removed 4 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 4 rows containing non-finite values (stat_smooth).

Warning: Removed 4 rows containing missing values (geom_point).

Warning: Removed 4 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 6 rows containing non-finite values (stat_smooth).
Warning: Removed 6 rows containing missing values (geom_point).

Warning: Removed 6 rows containing missing values (geom_point).

Warning: Ignoring unknown parameters: w, h
`geom_smooth()` using formula 'y ~ x'
Warning: Removed 10 rows containing non-finite values (stat_smooth).
Warning: Removed 10 rows containing missing values (geom_point).

Warning: Removed 10 rows containing missing values (geom_point).

Create long data frame

Disorders per Gender

Warning: Ignoring unknown parameters: binwidth

Warning: Ignoring unknown parameters: binwidth

100% bar per diagnosis with percentages of contribution per sample

Disorders per sample

Disorder_inc Sample n Prop
Anxiety disorders EDGI 38 0.00
Anxiety disorders GLAD 10098 0.66
Anxiety disorders NBR 1765 0.12
Anxiety disorders RAMP 3350 0.22
Depressive disorders EDGI 47 0.00
Depressive disorders GLAD 11455 0.63
Depressive disorders NBR 2756 0.15
Depressive disorders RAMP 3826 0.21
Eating disorders EDGI 61 0.03
Eating disorders GLAD 1384 0.65
Eating disorders NBR 195 0.09
Eating disorders RAMP 482 0.23
OCDs EDGI 18 0.01
OCDs GLAD 1627 0.68
OCDs NBR 145 0.06
OCDs RAMP 607 0.25
Psychotic disorders GLAD 401 0.68
Psychotic disorders NBR 37 0.06
Psychotic disorders RAMP 152 0.26
Bipolar disorders EDGI 3 0.00
Bipolar disorders GLAD 874 0.74
Bipolar disorders NBR 72 0.06
Bipolar disorders RAMP 240 0.20
PTSD EDGI 6 0.00
PTSD GLAD 1851 0.66
PTSD NBR 231 0.08
PTSD RAMP 701 0.25
Autism EDGI 6 0.01
Autism GLAD 384 0.61
Autism NBR 46 0.07
Autism RAMP 191 0.30
ADHD EDGI 2 0.01
ADHD GLAD 201 0.58
ADHD NBR 32 0.09
ADHD RAMP 114 0.33
Personality disorders EDGI 6 0.00
Personality disorders GLAD 951 0.72
Personality disorders NBR 76 0.06
Personality disorders RAMP 295 0.22
Controls EDGI 1 0.00
Controls GLAD 139 0.01
Controls NBR 9247 0.78
Controls RAMP 2538 0.21
Warning: Ignoring unknown parameters: binwidth

Disorder per ethnicity

Warning: Ignoring unknown parameters: binwidth

Age categories

Age bar plot for whole sample (not grouped) Summary dataset for age and display it as a bar plot

age_category n Prop order
16-18 1055 0.03 1
19-25 2617 0.08 2
26-35 5112 0.15 3
36-45 4781 0.14 4
46-55 6988 0.20 5
56-65 8015 0.23 6
66-70 3189 0.09 7
71-75 2010 0.06 8
76-80 544 0.02 9
81-85 116 0.00 10
86-90 27 0.00 11
91-100 6 0.00 12
100+ 5 0.00 13

Education

Pie charts

summary dataset with frequencies

highest_education n freq
GCSE/CSE 4224 0.12
NVQ 2148 0.06
A-levels 6742 0.20
University 19948 0.58

pie chart labels

Warning: Ignoring unknown parameters: face

Warning: Ignoring unknown parameters: face

Warning: Ignoring unknown aesthetics: face

pie chart - no labels

Ethnicity_collapsed

Pie charts

Ethnicity_collapsed n freq
European 30201 0.96
Mixed or multiple ethnic origins 571 0.02
Asian or Asian British 394 0.01
African or African British 159 0.01
Other 193 0.01

pie chart with only European label

Warning: Ignoring unknown parameters: face

Pie chart - no labels

Gender collapsed

Pie charts

Gender_collapsed n freq
Male 9698 0.28
Female 24387 0.71
Non-binary/Self-defined 380 0.01
Warning: Ignoring unknown parameters: face

Pie chart - only labelling two largest groups (male and female)

Warning: Ignoring unknown parameters: face

Pie chart - no labels

Compare the GLAD/EDGI and COPING subcohorts: EDGI, GLAD, NBR, RAMP

+++CH: Diagnosis, ethnicity, gender, age

age_category

Proportion tables

age_category Sample n Prop
16-18 EDGI 7 0.00
16-18 GLAD 341 0.01
16-18 NBR 10 0.00
16-18 RAMP 697 0.02
19-25 EDGI 17 0.00
19-25 GLAD 1537 0.04
19-25 NBR 291 0.01
19-25 RAMP 772 0.02
26-35 EDGI 21 0.00
26-35 GLAD 2850 0.08
26-35 NBR 1165 0.03
26-35 RAMP 1076 0.03
36-45 EDGI 8 0.00
36-45 GLAD 2459 0.07
36-45 NBR 1440 0.04
36-45 RAMP 874 0.03
46-55 EDGI 4 0.00
46-55 GLAD 2821 0.08
46-55 NBR 2644 0.08
46-55 RAMP 1519 0.04
56-65 EDGI 5 0.00
56-65 GLAD 1908 0.06
56-65 NBR 3901 0.11
56-65 RAMP 2201 0.06
66-70 EDGI 3 0.00
66-70 GLAD 433 0.01
66-70 NBR 1967 0.06
66-70 RAMP 786 0.02
71-75 GLAD 231 0.01
71-75 NBR 1265 0.04
71-75 RAMP 514 0.01
76-80 GLAD 59 0.00
76-80 NBR 334 0.01
76-80 RAMP 151 0.00
81-85 GLAD 13 0.00
81-85 NBR 61 0.00
81-85 RAMP 42 0.00
86-90 GLAD 1 0.00
86-90 NBR 13 0.00
86-90 RAMP 13 0.00
91-100 NBR 2 0.00
91-100 RAMP 4 0.00
100+ NBR 3 0.00
100+ RAMP 2 0.00
Sample age_category n Prop
EDGI 16-18 7 0.11
EDGI 19-25 17 0.26
EDGI 26-35 21 0.32
EDGI 36-45 8 0.12
EDGI 46-55 4 0.06
EDGI 56-65 5 0.08
EDGI 66-70 3 0.05
GLAD 16-18 341 0.03
GLAD 19-25 1537 0.12
GLAD 26-35 2850 0.23
GLAD 36-45 2459 0.19
GLAD 46-55 2821 0.22
GLAD 56-65 1908 0.15
GLAD 66-70 433 0.03
GLAD 71-75 231 0.02
GLAD 76-80 59 0.00
GLAD 81-85 13 0.00
GLAD 86-90 1 0.00
NBR 16-18 10 0.00
NBR 19-25 291 0.02
NBR 26-35 1165 0.09
NBR 36-45 1440 0.11
NBR 46-55 2644 0.20
NBR 56-65 3901 0.30
NBR 66-70 1967 0.15
NBR 71-75 1265 0.10
NBR 76-80 334 0.03
NBR 81-85 61 0.00
NBR 86-90 13 0.00
NBR 91-100 2 0.00
NBR 100+ 3 0.00
RAMP 16-18 697 0.08
RAMP 19-25 772 0.09
RAMP 26-35 1076 0.12
RAMP 36-45 874 0.10
RAMP 46-55 1519 0.18
RAMP 56-65 2201 0.25
RAMP 66-70 786 0.09
RAMP 71-75 514 0.06
RAMP 76-80 151 0.02
RAMP 81-85 42 0.00
RAMP 86-90 13 0.00
RAMP 91-100 4 0.00
RAMP 100+ 2 0.00

Cross tables

age_category n_EDGI n_GLAD n_NBR n_RAMP Prop_EDGI Prop_GLAD Prop_NBR Prop_RAMP
16-18 7 341 10 697 0.11 0.03 0.00 0.08
19-25 17 1537 291 772 0.26 0.12 0.02 0.09
26-35 21 2850 1165 1076 0.32 0.23 0.09 0.12
36-45 8 2459 1440 874 0.12 0.19 0.11 0.10
46-55 4 2821 2644 1519 0.06 0.22 0.20 0.18
56-65 5 1908 3901 2201 0.08 0.15 0.30 0.25
66-70 3 433 1967 786 0.05 0.03 0.15 0.09
71-75 NA 231 1265 514 NA 0.02 0.10 0.06
76-80 NA 59 334 151 NA 0.00 0.03 0.02
81-85 NA 13 61 42 NA 0.00 0.00 0.00
86-90 NA 1 13 13 NA 0.00 0.00 0.00
91-100 NA NA 2 4 NA NA 0.00 0.00
100+ NA NA 3 2 NA NA 0.00 0.00
age_category EDGI GLAD NBR RAMP
16-18 7 341 10 697
19-25 17 1537 291 772
26-35 21 2850 1165 1076
36-45 8 2459 1440 874
46-55 4 2821 2644 1519
56-65 5 1908 3901 2201
66-70 3 433 1967 786
71-75 0 231 1265 514
76-80 0 59 334 151
81-85 0 13 61 42
86-90 0 1 13 13
91-100 0 0 2 4
100+ 0 0 3 2

Chi square test

Warning in chisq.test(.): Chi-squared approximation may be incorrect
statistic p.value parameter method
5856.536 0 36 Pearson’s Chi-squared test
Warning: Ignoring unknown parameters: binwidth

Warning: Ignoring unknown parameters: binwidth

Ethnicity_collapsed

Proportion tables

Ethnicity_collapsed Sample n Prop
European EDGI 60 0.00
European GLAD 12090 0.35
European NBR 9927 0.29
European RAMP 8124 0.24
Mixed or multiple ethnic origins EDGI 3 0.00
Mixed or multiple ethnic origins GLAD 276 0.01
Mixed or multiple ethnic origins NBR 118 0.00
Mixed or multiple ethnic origins RAMP 174 0.01
Asian or Asian British EDGI 1 0.00
Asian or Asian British GLAD 98 0.00
Asian or Asian British NBR 126 0.00
Asian or Asian British RAMP 169 0.00
African or African British GLAD 37 0.00
African or African British NBR 46 0.00
African or African British RAMP 76 0.00
Other EDGI 1 0.00
Other GLAD 122 0.00
Other RAMP 70 0.00
NA GLAD 30 0.00
NA NBR 2879 0.08
NA RAMP 38 0.00
Sample Ethnicity_collapsed n Prop
EDGI European 60 0.92
EDGI Mixed or multiple ethnic origins 3 0.05
EDGI Asian or Asian British 1 0.02
EDGI Other 1 0.02
GLAD European 12090 0.96
GLAD Mixed or multiple ethnic origins 276 0.02
GLAD Asian or Asian British 98 0.01
GLAD African or African British 37 0.00
GLAD Other 122 0.01
GLAD NA 30 0.00
NBR European 9927 0.76
NBR Mixed or multiple ethnic origins 118 0.01
NBR Asian or Asian British 126 0.01
NBR African or African British 46 0.00
NBR NA 2879 0.22
RAMP European 8124 0.94
RAMP Mixed or multiple ethnic origins 174 0.02
RAMP Asian or Asian British 169 0.02
RAMP African or African British 76 0.01
RAMP Other 70 0.01
RAMP NA 38 0.00

Cross tables

Ethnicity_collapsed n_EDGI n_GLAD n_NBR n_RAMP Prop_EDGI Prop_GLAD Prop_NBR Prop_RAMP
European 60 12090 9927 8124 0.92 0.96 0.76 0.94
Mixed or multiple ethnic origins 3 276 118 174 0.05 0.02 0.01 0.02
Asian or Asian British 1 98 126 169 0.02 0.01 0.01 0.02
Other 1 122 NA 70 0.02 0.01 NA 0.01
African or African British NA 37 46 76 NA 0.00 0.00 0.01
NA NA 30 2879 38 NA 0.00 0.22 0.00
Ethnicity_collapsed EDGI GLAD NBR RAMP
European 60 12090 9927 8124
Mixed or multiple ethnic origins 3 276 118 174
Asian or Asian British 1 98 126 169
African or African British 0 37 46 76
Other 1 122 0 70
NA 0 30 2879 38

Chi square test

Warning in chisq.test(.): Chi-squared approximation may be incorrect
statistic p.value parameter method
5106.32 0 15 Pearson’s Chi-squared test
Warning: Ignoring unknown parameters: binwidth

Warning: Ignoring unknown parameters: binwidth

Gender_collapsed

Proportion tables

Gender_collapsed Sample n Prop
Male EDGI 1 0.00
Male GLAD 2279 0.07
Male NBR 5641 0.16
Male RAMP 1777 0.05
Female EDGI 62 0.00
Female GLAD 10127 0.29
Female NBR 7416 0.22
Female RAMP 6782 0.20
Non-binary/Self-defined EDGI 2 0.00
Non-binary/Self-defined GLAD 247 0.01
Non-binary/Self-defined NBR 39 0.00
Non-binary/Self-defined RAMP 92 0.00
Sample Gender_collapsed n Prop
EDGI Male 1 0.02
EDGI Female 62 0.95
EDGI Non-binary/Self-defined 2 0.03
GLAD Male 2279 0.18
GLAD Female 10127 0.80
GLAD Non-binary/Self-defined 247 0.02
NBR Male 5641 0.43
NBR Female 7416 0.57
NBR Non-binary/Self-defined 39 0.00
RAMP Male 1777 0.21
RAMP Female 6782 0.78
RAMP Non-binary/Self-defined 92 0.01

Cross tables

Gender_collapsed n_EDGI n_GLAD n_NBR n_RAMP Prop_EDGI Prop_GLAD Prop_NBR Prop_RAMP
Male 1 2279 5641 1777 0.02 0.18 0.43 0.21
Female 62 10127 7416 6782 0.95 0.80 0.57 0.78
Non-binary/Self-defined 2 247 39 92 0.03 0.02 0.00 0.01
Gender_collapsed EDGI GLAD NBR RAMP
Male 1 2279 5641 1777
Female 62 10127 7416 6782
Non-binary/Self-defined 2 247 39 92

Chi square test

Warning in chisq.test(.): Chi-squared approximation may be incorrect
statistic p.value parameter method
2455.701 0 6 Pearson’s Chi-squared test
Warning: Ignoring unknown parameters: binwidth

Warning: Ignoring unknown parameters: binwidth

highest_education

Proportion tables

highest_education Sample n Prop
GCSE/CSE EDGI 7 0.00
GCSE/CSE GLAD 1078 0.03
GCSE/CSE NBR 1876 0.06
GCSE/CSE RAMP 1263 0.04
NVQ EDGI 3 0.00
NVQ GLAD 918 0.03
NVQ NBR 1227 0.04
A-levels EDGI 20 0.00
A-levels GLAD 2713 0.08
A-levels NBR 2320 0.07
A-levels RAMP 1689 0.05
University EDGI 35 0.00
University GLAD 7515 0.23
University NBR 6983 0.21
University RAMP 5415 0.16
Sample highest_education n Prop
EDGI GCSE/CSE 7 0.11
EDGI NVQ 3 0.05
EDGI A-levels 20 0.31
EDGI University 35 0.54
GLAD GCSE/CSE 1078 0.09
GLAD NVQ 918 0.08
GLAD A-levels 2713 0.22
GLAD University 7515 0.61
NBR GCSE/CSE 1876 0.15
NBR NVQ 1227 0.10
NBR A-levels 2320 0.19
NBR University 6983 0.56
RAMP GCSE/CSE 1263 0.15
RAMP A-levels 1689 0.20
RAMP University 5415 0.65

Cross tables

highest_education n_EDGI n_GLAD n_NBR n_RAMP Prop_EDGI Prop_GLAD Prop_NBR Prop_RAMP
GCSE/CSE 7 1078 1876 1263 0.11 0.09 0.15 0.15
NVQ 3 918 1227 NA 0.05 0.08 0.10 NA
A-levels 20 2713 2320 1689 0.31 0.22 0.19 0.20
University 35 7515 6983 5415 0.54 0.61 0.56 0.65
highest_education EDGI GLAD NBR RAMP
GCSE/CSE 7 1078 1876 1263
NVQ 3 918 1227 0
A-levels 20 2713 2320 1689
University 35 7515 6983 5415

Chi square test

Warning in chisq.test(.): Chi-squared approximation may be incorrect
statistic p.value parameter method
1207.855 0 12 Pearson’s Chi-squared test

Bar plots

Warning: Ignoring unknown parameters: binwidth

Warning: Ignoring unknown parameters: binwidth

key_worker

Proportion tables

key_worker Sample n Prop
No key worker EDGI 42 0.00
No key worker GLAD 7249 0.21
No key worker NBR 8660 0.25
No key worker RAMP 5841 0.17
Key worker EDGI 23 0.00
Key worker GLAD 5366 0.16
Key worker NBR 4414 0.13
Key worker RAMP 2765 0.08
NA GLAD 38 0.00
NA NBR 22 0.00
NA RAMP 45 0.00
Sample key_worker n Prop
EDGI No key worker 42 0.65
EDGI Key worker 23 0.35
GLAD No key worker 7249 0.57
GLAD Key worker 5366 0.42
GLAD NA 38 0.00
NBR No key worker 8660 0.66
NBR Key worker 4414 0.34
NBR NA 22 0.00
RAMP No key worker 5841 0.68
RAMP Key worker 2765 0.32
RAMP NA 45 0.01

Cross tables

key_worker n_EDGI n_GLAD n_NBR n_RAMP Prop_EDGI Prop_GLAD Prop_NBR Prop_RAMP
No key worker 42 7249 8660 5841 0.65 0.57 0.66 0.68
Key worker 23 5366 4414 2765 0.35 0.42 0.34 0.32
NA NA 38 22 45 NA 0.00 0.00 0.01
key_worker EDGI GLAD NBR RAMP
No key worker 42 7249 8660 5841
Key worker 23 5366 4414 2765
NA 0 38 22 45

Chi square test

Warning in chisq.test(.): Chi-squared approximation may be incorrect
statistic p.value parameter method
332.6547 0 6 Pearson’s Chi-squared test

Bar plots

Warning: Ignoring unknown parameters: binwidth

Warning: Ignoring unknown parameters: binwidth